arXiv: 1502.02564v2 [hep-th] 25 Aug 2015 


Note on the Gauss-Bonnet braneworld scenario 


D. Bazeia^, A. Lobao^, L. Losano^, R. Menezes^’"^, A. Yu. Petrov^ 

^ Departamento de Fisica, Universidade Federal da Paraiha, 

58051-970, Jodo Pessoa, Paraiba, Brazil 
‘^Eseola Teeniea de Saude de Cajazeiras, 

Universidade Federal de Campina Grande, 58900-000 Cajazeiras, PB, Brazil 
^Departamento de Cieneias Exatas, Universidade Federal da Paraiba, 

58297-000 Rio Tinto, PB, Brazil and 
'^Departamento de Fisiea, Universidade Federal de Campina Grande, 

58109-970 Campina Grande, PB, Brazil 
(Dated: August 26, 2015) 

Abstract 

In this work we deal with the presence of braneworld solutions in a five-dimensional space-time 
with a single extra spatial dimension of infinite extent. The braneworld scenario is built under the 
presence of a single real scalar field, and we modify the gravity sector to include generic function 
of the Gauss-Bonnet term. We study several specific models, and we construct exact braneworld 
solutions, in particular for including the Gauss-Bonnet term at first and second order power. As 
an interesting result, we show that the brane tends to split for a specific modification in the gravity 
sector, in the presence of non constant Gauss-Bonnet term. 

PACS numbers: 04.50.-h, ll.25.-w 
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The study of gravity is certainly one of the most fundamental research lines in the modern 
theoretical high energy physics. The problem of explaining the cosmic acceleration discov¬ 
ered in [T], together with search for a consistent quantum gravity model, called a great deal 
of attention to the study of modified theories of gravity. The crucial role was played by the 
paper [2] where the idea of noncompact extra dimensions, which gave rise to wide application 
of the brane concepts in gravity and cosmology, has been introduced. 

Following the brane concept, our four-dimensional space-time is treated as a brane em¬ 
bedded to some higher-dimensional space-time, that is, the bulk where the gravitational field 
can propagate, whereas the fields of fundamental interactions are confined to the brane [3] . 
In the standard scenario, the higher-dimensional gravity was supposed to be described by 
the usual Einstein-Hilbert action. However, a very natural development of the idea consisted 
in modifying of the gravity Lagrangian implemented through the replacement of the scalar 
curvature R by some other scalar depending on the gravitational helds and involving its 
higher derivatives. Within the quantum held theory, the interest to such a modihcation of 
the gravity has been arisen by the seminal paper |1] where the higher-derivative extensions 
of gravity were considered as a way to achieve a renormalizable gravity theory. Within the 
classical gravity, study of the f{R) models was developed as an aim for a possible explana¬ 
tion of a cosmic acceleration |5]. Many issues related to the different aspects of the f{R) 
gravity, especially the exact solutions in this theory, have been studied in [B]. Within the 
brane context, the f{R) gravity has been applied for the hrst time in [7]. 

However, the f{R) gravity does not include all possibilities for extensions of gravity. One 
of the most studied extensions is the Lovelock gravity, such that the gravity action is now 
represented as a series, with its lowest (zero) order being a cosmological term, the first order 
is the scalar curvature, and the second one is the Gauss-Bonnet term [8]. As a result, the 
Gauss-Bonnet gravity naturally emerges, in which the gravity Lagrangian is a sum of the 
usual Einstein-Hilbert Lagrangian with the Gauss-Bonnet scalar or its function. It is well- 
known that in the four-dimensional space-time, the Gauss-Bonnet term is topological. Some 
aspects of the Gauss-Bonnet gravity have been considered earlier in [9]. In m, this theory 
has been applied within the cosmic acceleration context. It should be mentioned, that the 
higher-order contractions of the Riemann curvature tensor, and hence higher orders in the 
Gauss-Bonnet term, naturally arise in the low-energy limit of the string theory |T0] . 

Within this paper, we are going to study the Gauss-Bonnet gravity within the brane con¬ 
text. Some earlier studies in this direction have been performed in [T^ 113] . where, however, 
only the simplest version of the Gauss-Bonnet gravity has been considered, that is, only 
the Gauss-Bonnet term G itself is added to the gravity action. At the same time, actually 
more sophisticated manners to include the Gauss-Bonnet term are carried out, implying 
thus in the theories of f{G) gravity and f{R,G) gravity which are intensively discussed 
within the cosmic acceleration problem, for a review see f.e. [H] and references therein. 
It is worth mentioning that f{G) gravity has been tested within the cosmological observa¬ 
tions in na, where the restrictions on the numerical parameters of the action have been 
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obtained. Besides of the cosmological studies, the stability issues also have been discussed 
for the Gauss-Bonnet gravity [1^, where, moreover, it was shown that just for the case 
/(G) = G, the Gauss-Bonnet gravity displays the instabilities, which certainly implies a 
suggestion that namely a theory described by a more sophisticated function /(G) can be 
free of instabilities. Therefore, it seems natural to apply these more generalized models also 
within the braneworld context. 

Our aim will consist principally in applying the hrst-order formalism which has manifested 
itself as a very useful tool for solving the nonlinear equations of motion, especially modihed 
Einstein equations [7j. Using this method, we will hnd some new exact solution for the 
Gauss-Bonnet brane. The exact solutions are then used to investigate the prohle of the 
energy density, searching in particular for the brane splitting feature, as it has appeared 
before in dUEH] and more recently in |T9]. In the recent work [19], one noticed that the 
splitting appears due to modihcation in the geometry, with the inclusion of the B? term 
in the Einstein-Hilbert action. In the current work we modify the geometry including the 
Gauss-Bonnet term, so we believe that the splitting is not excluded, being an effect unseen 
in previous investigations. 

We start with the following action describing the Gauss-Bonnet brane (cf. |13jh 

S = j d^xdy^g Qi? + ^/(G) + . (1) 

Here y is the extra coordinate, Cg is the source Lagrange density. As usual, G stands for 
the Gauss-Bonnet term, which is given by 

G = R^- + R^,xpR>^''^P. (2) 

In the present work, we take the line element in the form 

ds^ = e^'^^^^rjabdx°'dx^ — dy"^, (3) 

with A representing the warp function. This stands for AdS^ geometry, and here we consider 
the source as a scalar held, specihed by 

c. = \g,.d^‘4>d"4> - V{4>). (4) 

The modihed Einstein equation looks like [13] 

G^, = 2T,,+ ^-g^J{G)-2F{G)RR,, + 4F{G)R^^R,x-2F{G)R^xpaR,"^^ 

- AF{G)R^p^,RP^ + 2i?V^V.F(G) - 2Rg^,V‘^F{G) - 4i?^V.VpF(G) 

- 4i?eVpVpF(G) + 4R,,V^F{G) + 4g,,R^^VxVpF{G) - AR^^xpV^V^FiG) 

= (5) 

Here we are using F{G) = T)(i/ is the energy-momentum tensor of the source held, 

and Gpy = R^y — \Rg^v is the usual Einstein tensor. Since G is a scalar depending only 
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on y (= 0 : 4 ), by symmetry the last term AR^yXpV^V^FiG) = 0 ; unfortunately, however, this 
does not simplify the expression too much. For the sake of simplicity we choose in this 
paper f{G) = aG”, with n > 1. We use Latin indices a,b,c,d = 0...3, and Greek ones 
/iz/, Ap = 0...4. Also, A' = It is clear that for f{G) = 0 the usual equations of motion 
[Is] are restored. 

Direct calculations show that for the braneworld metric ([^ one has 


Rabcd = e^^A.'‘^{r]acr]bd - VadVbc)] RibAd = -e^^(A" + 

Rbd = e^\d(A" + 4A'2); i?44 = -4 (A" + A'2); R = 8 A" + 20 A'^; 

Rp^Rf^^ = 80 A'^ + 20A"2 + 64A"A'2; Rp^xpR^'^^^ = 40 A'^ + 16 A"^ + 32 A" A'^; 

G = 24[4A"A'2+ 5A'^]; 

Tab = Vabil^P'^ + D(0))e2^; T 44 = ^0'^ - D(0). ( 6 ) 

For the usual Einstein tensor Gpi, one hnds in this metric 


Gab — —S{A" + 2A'‘^)rjab&^^] G 44 — 6 A^^; GaA — 0. (7) 


Let us now look for some possibilities to solve the equations ([^. The most complicated 
terms are the derivatives of /(G). As a hrst step, one can verify that the expression for the 
Gauss-Bonnet term is correct. Indeed, let us suggest that the metric in the D-dimensional 
bulk space-time looks like ([^, but with a, 6 = 0...D — 2. Repeating all calculations, one gets 

G = (D - 1)(D - 2)(D - 3)[DA'^ + AA" ( 8 ) 


At D = 5, one indeed yields G = 120A''^-|-96A"A'^. At D = 4, one has G 4 = 24A'‘^-|-24A"A'^; 
in this case, the contribution to the action, that is, ^/\^\G 4 ^ = 24e^^(A'^-f A"A'^) = 8 (e^^A'^)' 
(recall that \/|^ = is a total derivative, thus conhrming the known fact that the 

Gauss-Bonnet Lagrangian is a total derivative in four dimensions. Also, it vanishes in lower 
dimensions. 

The next step is to suggest that the Gauss-Bonnet term G is constant, i.e. 


G 


4A"A'2 + 5A'^ = bB. 


In this case, for the positive B = 6 ^, one gets 


y + G 


4 

5 


dA' A'2 
b^-A'^' 


and for the negative B = — 6 ^, it gives 


y + G 


4 r dA' A '2 

5 J 62 + A'^ ■ 


(9) 


( 10 ) 


( 11 ) 


4 




In principle, the explicit forms of A{y) can be found for both these factors, although it is 
ugly: for B = b'^, one has 


y + C 


■-( — 
5 V4y^ 


In 


A' + Vb 
A' -Vb 


1 A'\ 

—-= arctan —= , 

2^fb VbB 


( 12 ) 


and for B = —6^, one has an even more complicated expression. However, the result 
evidently exists. So, let us write the equations of motion for this case: 


G 




+ \9^.uf{G) - 2F{G)RR^, + 4F{G)R^^R,x - 2F(G')/2^Apa V 
4F{G)R^p^,R>^F (13) 


Here f{G) and F{G) are constants. So, we have for the cases pz/ = ab and yu = 44, 
respectively: 

-3(H" + 2H'2) = (0'2 + 2H(0)) + ^/(G)-F(G)(48H'^ + 36H"H'2); (14) 

6H'2 = _ 2V{(j)) - i/(G) + F(G')(48H'^ + 48H"H'2), 

or, as is the same after replacement of the hrst equation by the sum of the two equations. 


-3H" = 20'2 + 12F(G)H"H'2; 

- 21/(0) - ^f{G) + A8F{G){A'^ + A"A'^). (15) 

So, we can use some choice for the warp factor consistent with the constant G and arrive 
at the system on 0 and 1/(0). We note that within this methodology, the potential 1/(0), 
instead of being introduced as a basic characteristic of the theory, turns out to be a variable 
to be determined. However, the assumptions of this method, imposing restrictions on the 
potential mm, are necessary to employ the power of the first-order formalism allowing to 
obtain exact solutions for a certain class of potentials while, in a generic case, the solutions 
can only be found numerically. Therefore, within this approach, we solve a kind of inverse 
problem, allowing to hnd the potential for the known warp factor. Moreover, we demonstrate 
how the explicit form of the potential can be found. 

Let us substitute G = 120/1, with D being a constant. Also, we suppose that if A' = 
A'{y), one can use y = y{A') as well, therefore one has 


0 = 5 


D-A'^)\ 00 
4A'2 J dW 


( 16 ) 


We substitute A' = 2 ;. Thus, our system ([15|) implies 


- 15 


D-z^ 




D-z^ 

4:Z‘^ 


4\ 2 


+ mF{D){D - z'^); 


6z^ = 25 


D-z^ 

4^2 


^ ) - 21/(0) - ^f{D) + 48F{D){5D - z^) ,(17) 
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where F{D) = F{G{D)) = F{120D), similarly, f{D) = f{120D). Resolving these equa¬ 
tions, we hnd 


/^ 
\dz J 

v(zm 


1 ^ ~ 4^2 

g0(l + 4W)^; 

§ Y(1 + 44^(45)*-^ - - \hD) + 

48F(D)(f}D - z‘) . 


(18) 


So, in principle, the held and the potential are found in terms of the auxiliary variable 
At the same time, we can abandon the restriction of the constant G. In this case, the 
Hab and Fl/^^ terms take the form 

Hab = Vabe^^ “ 12F{G)A'^ {4A'^ + 3A") + 24F'{G)A' {A'^ + A') 

+ AF'\G) (2A'2 - A') 1; 

Hm = -\fiG) + A8F{G)A^ {a" + A^) + 1QF'{G)A {A" - 2A‘^) . (19) 

We consider these equations for the particular case f{G) = aG'^ which yields 

Hab = Vabe^^[^aG^ - 12anA^ {4A^ + 3A") G^-^ 

+Aan{n - 1) \qG'A {A^ + A") + {2A'^ - A") G"] G'”-^ 

+4an{n - l)(n - 2)G'^ {2A^ - A") 

//44 = + 48 anA '2 (W' + A^) + 16an(n - 1)G'A' {A' - 2A^) G”-^. (20) 

Summing up the equations for T (the T is an object dehned from the relation Tab = Vab^^^T) 
and T 44 , we hnd 

V" = -^A" - - li/44. (21) 

Therefore, one has 

= --A' - QanA^A'G^-^ 

-2an{n - 1 ) [ {2A‘^ - A") G" - 2A (A'^ - 5A") G'j G”-^ 

-2an{n - 1)(n - 2)G'‘^ (2A^ - A") G”-^ (22) 

It is instructive to give here the explicit expressions for the potential and energy density. 
To do it, we note that the modihed Einstein equations will look like 

- 3(A" + 2A^) = (0'2 + 2V) + if; 

6A'2 = -2R)+ii44, (23) 


6 






FIG. 1: The warp factor (28), depicted for B = 0.5 (red, solid line), B 
and B = 20.0 (blue, dotted-dashed line). 


3.0, (brown, dashed line) 


where H can be read off from Eq. (20). Subtracting the second equation from the first one, 
we find the potential: 


^ (-3(dl" + AA'^)-H + Hu) , 

whose explicit form is 


(24) 


1 /( 0 ) = -^A" - 3 / 1'2 - + 3anA'^ {8A'^ + 7A") 

-an{n - 1) (/!' + 2A") G' + (2A'^ - A") G"Jg^-^ 

-an(n - l)(n - 2)G'^ (2/1'^ - A") G^-\ (25) 


The energy density in our case is given by the Tqq component of the energy-momentum 
tensor of the matter 

p = (7r + r(^)) . (26) 

Let us now consider some specihc cases. 

a. First example. For n = 1 and a = —1/4, we hnd that 

= -^A"+^A'^A” . (27) 


We choose a standard ansatz (which at large \y\ tends to A{y) = —B\y\): 


A{y) = 51n[sech(?/)], 


(28) 


where B > 0. We plot the warp function for the above warp function in Fig. We 

note that the warp factor becomes more and more localized around y = 0, as B increases to 
larger and larger values. 


We use the Eq. (27) to get 


l/2 


= ^ (l - B'^) sech^{y) + ^sech^(|/). 


(29) 
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Moreover, in the case B = 1, we find 


t>' = ±\l'^sech.‘^{y). 


The solution of this equation is kink-like: 


(30) 


</>(?/) 


± 



tanh(j/). 


(31) 


b. Second example. 
now we consider n = 2. 


We use the same ansatz for the warp function, as in Eq. (28), but 
We get 


_3_ 4c, (sA'^A'^G + (2A'2 - A") G" - 2^4' - 5A") G") . (32) 

Using the expression for the Gauss-Bonnet term (|^, we arrive at 

= --A” + 96a(25A'^A” -32A'^A” -SA^A”” + SA”^ -316A^A”^ 

-36A'^A"^ - 68A^A"A” + AA^A”A”” + 2AAA”^A”^ . (33) 

Using the warp function ( [2^ , we find 

{1 + 64aB\32 + 64B - 25B‘^)) sech\y) 

+96aB^ (16 - 248S - 5085^ + 75B^) sedi^{y) 

-9QaB^ (88 - 4365 - 8245^ + 755^) sech®(|/) 

+480a5^ (16 - 445 - 765^ + 55^) sech®(?/), (34) 


Now, we choose 


a = 


which allows to write Eq. (34) as 


1 

645^(32 + 645 - 2552) ’ 


cj)''^ = Oi sech'^(j/) — 02 sech®(j/) -|- 03 sech®(j/), 


where 


Oi 


02 


03 


Let us rewrite Eq. (36) as 


3(16 - 2485 - 50852 ^ 75 ^ 3 ^ 
2(32 + 645 - 2552 ) 

3(88 - 4365 - 8245^ + 755^) 
2(32 + 645 - 2552 ) 
15(16 - 445 - 7652 ^ 5 ^ 3 ^ 
2(32 + 645 - 2552 ) ■ 


0 '^ = 03 sech^(?/) ^ sech 2 (|/) -f sech'^(?/) 

L03 03 


(35) 

(36) 

(37) 

(38) 

(39) 


( 40 ) 












FIG. 2: The solutions (46) (red, solid line), l/10x([47|) (brown, dashed line), and l/2x(48) (blue, 
dotted-dashed line). 


Now, we define the constants 


with the restriction 


Oi 2 

— = r 

Os 


^ = 2r, 

03 


± / — - 


02 


O3 2o3 


that is, 


± 


'I 6 - 2485 - 50852 ^ 7553 


4365 - 82452 ^ 75^3 


(41) 


(42) 


(43) 


5(16 - 445 - 7652 + 553 ) io(i 6 - 445 - 7652 + 553 ) • 

One can show that the possible solutions correspond only to the negative sign in this 
equation which possesses six solutions, with two of them being essentially complex. From 
the remaining four solutions, only three satisfy 0'^ > 0, they are: 5i = 0.4524, B 2 = 2.9561 


and 53 = 20.0107. Therefore, if Eq. (43) is satisfied, from Eq. (40) we have 


(j)' = ifycfy secfy(|/) r — sech^(|/) 


which has the general solution 


1 


(t){y) = 3r - 2 - sech^{y) tanh(|/). 

From that, for permitted values of 5, we obtain 

My) = ± (0.5325tanfy(?/)- 3.8757tanh( 2 /)) for 
My) = ± (14.0826 tanfy(j/) — 84.3242 tanh(j/)) for 
My) = ± (0.9165tanfy(l/) — 10.9778 tanh(|/)) for 


5 = 5i, 
5 = 52 , 
5 = 53 , 


(44) 


(45) 


(46) 

(47) 

(48) 


which are depicted in Fig. (|^. We note that although the solutions vary similarly, the 
amplitude of 02 (l/) is much greater then it is in the other two cases. 

In Fig. 1^ we depict the energy density obtained from (26), corresponding to each one 
of the above solutions. We see that the brane has the standard behavior for 5 = 53 , but 
it splits significantly for B = B 2 . The splitting behavior of the brane was studied before 
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FIG. 3: The energy densities for the solutions (46) (red solid line), (47)(brown dashed line), and 
(48) (blue dotted-dashed line). 


in m [18] . as an effect driven by the source Lagrange density. It was also identified in 
[in], through modihcation of the geometry, shown to appear from the inclusion of the 
term, in the case of non constant curvature. Here it also appear through modification of the 
geometry, with the inclusion of the term in the Einstein-Hilbert action, for non constant 
Gauss-Bonnet term. Although the splitting is similar to the case found in Ref. [19], in 
the present case it is much more evident, making the core of the brane to behave as it is 
asymptotically. 

In these cases, the potential can be found explicitly. For the case n = 1 it is given by 

^(0(l/)) = + |Sech2(|/) + |sech^(|/). (49) 

The solution for the field is 


0 = 



tanh(?/) ^ y 


arctanh 



Therefore, we can explicitly express the potential V (0) as a function of the field: 

For the case n = 2 the potential looks like 

V{(j){y))=—3 + 2160a + (^2 11664ajsech^(|/) — 62736asech"^(|/) -|- 

+77808asech®(|/) - 28512a sech®(|/). 


(50) 


(51) 


(52) 


where a is given by Eq. (35). In this case we have 


1 /( 0 ) = arctanh(s(0)), 


(53) 


where s(0) satishes the algebraic equation 

as^(0) — 6s(0) — 0 = 0 , 


(54) 


with a and b are the numerical parameters which can be read off from Eqs. ( [^47[48 ). 
Therefore, one has the potential 

1/(0) = ^ -|- 384a — ^5568a -|- ^ js^(0) — 384a s^(0) + 36240a s®(0) — 28512a s®(0), (55) 
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so we conclude that we succeeded to obtain the potential in the form V = V{(j)). 

We can go on and use in (22) the ansatz A' = 0' = as in [7j. We take 

into account that the equation of motion of the scalar held is the same as in the usual f{R) 
brane case [7] 


whose solution is 




VU) = — (mD^ + 
2m 


(56) 


(57) 


just as in [7]; so, we hnd the exponential potential. In this case, the scalar held yields the 
Liouville-like form: 

(j)=-—\n{mD{y-yo)), =-^0. (58) 

m D 


These ansatzes allow to reduce the equation (22) to a purely algebraic, although rather 
complicated, equation on the parameter m which can be solved. Its explicit form is 


{D^ + -mCD)e^"''l’ = -e'‘”’"*[24C“(5C + 4mB)]"-'x 


(59) 


X a[6nmDC^ + 2n{n — l)[22m^C^D^ — Gm^CD^ — 2mC^D]] 


We see that here one meets two situations: hrst, one can have n = 1/2 (which is rather 
exotic - it means that the Gauss-Bonnet term enters the action with a fractional degree), 
second, the r.h.s. and the l.h.s. of this equation are separately equal to zero. In both 
cases one rests with some algebraic equations for coefficients C and D, arriving at the hnal 
expression D = —3mC/2, with 97m^(n — 1) = 6; thus, for the given n, we get a fractional 
m. 

In summary, in this work we considered the Gauss-Bonnet braneworld model. For this 
theory, we used the hrst-order formalism which allowed us to hnd some exact solutions. 
First, we succeeded to hnd the scalar held and its potential in an implicit form, in terms 
of some auxiliary variable, for the case of the constant Gauss-Bonnet term. Second, we 
abandoned the condition that the Gauss-Bonnet term is constant, and in this case we found 
the solutions for the specihc forms of the function f{G), that is, f{G) = aG"' with n = 1 and 
n = 2. Also, we found the solution corresponding to the exponential potential. All these 
solutions are not compatible with the constant scalar curvature, therefore, the anti-de-Sitter 
space is ruled out, and the supersymmetric extension of the theory seems to be impossible. 

As a particularly interesting behavior, we noted the splitting of the brane for non constant 
Gauss-Bonnet term in the case of n = 2, that is, for f{G) = aG^] see Fig. This behavior 
appears from modihcation of the geometry, and it is similar to the splitting found before in 
[in], but here it is much more evident for B = B 2 = 2.9561, with the core of the solution 
behaving as it is in the region far away from the brane. 


11 




An interesting study concerns stability of the braneworld scenarios that we investigated 
above. We can follow two distinct approaches, one investigating how the tensorial fluctua¬ 
tions in the metric behave, and the other, adding small parameters to control the behavior 
of the extra terms one has introduced in the current study. These and other similar issues 
are presently under consideration, and we hope to report on them in the near future. 
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